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^ ■ ABSTRACT 

I We report the discovery of a low-mass (0.26±0.02, 0.15±0.01 Mq) pre-main-sequence 

'— 'i eclipsing binary with a 5.3 day orbital period. JW 380 was detected as part of a 

high-cadence time-resolved photometric survey (the Monitor project) using the 2.5 m 
Isaac Newton Telescope and Wide Field Camera for a survey of a single field in the 
Orion Nebula Cluster (ONC) region in V and i bands. The star is assigned a 99 per 
' cent membership probability from proper motion measurements, and radial velocity 

, observations indicate a systemic velocity within la of that of the ONC. Modelling 

psj ■ of the combined light and radial velocity curves of the system gave stellar radii of 

, — ! I 1.19lQ;°g R© and 0.90lg Rq for the primary and secondary, with a significant third 

. light contribution which is also visible as a third peak in the cross-correlation functions 

' used to derive radial velocities. The masses and radii appear to be consistent with 

I stellar models for 2 — 3 Myr age from several authors, within the present observational 

• • ■ errors. These observations probe an important region of mass-radius parameter space, 

, ' where there are currently only a handful of known pre-main-sequence eclipsing binary 

, systems with precise measurements available in the literature. 

i Key words: stars: individual: JW 380 - open clusters and associations: individual: 

■ ■ ■ ' Orion Nebula Cluster - binaries: eclipsing - stars: pre-main-sequence - surveys. 
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1 INTRODUCTION 

Detached eclipsing binaries (EBs) provide one of the most 
precise ( ^ 2 per cent) and accurate (largely model and dis- 
tance independent) methods for measurement of fundamen- 
tal stellar properties (particularly, masses and radii). These 
can be used to place stringent constraints on stellar evolu- 
tion models. 

On the pre-main-sequence (PMS), such constraints are 
presently extremely scarce below 1 M0, and to our knowl- 
edge there are only 5 known PMS EB systems in this mass 
range. These are a 1.0, 7 Mp EB fStassu n et al] 120041') 
and a 1.27, 0.93 Mq EB l|Covino et al.,,2004i : iCovino et all 



both in Orion (the former is thought to belong to 
the Ori Ic association), with corresponding age ~ 5 — 
10 Myr, an M- dwarf EB in NGC 1647 (~ 150 Myr; 
iHebb et al.ll2006l '). and t wo EBs in the ONC: a brown dwarf- 
brown dwarf system (Stassun. Mathieu fc Valentil l2006l : 
IStassun. Mathieu fc Valentil 120071 ) and a 0.4. 0.4 Mr.^ M- 
dwarf EB (Cargile, Stassun & Mathieu, submitted). Com- 
parison of the NGC 1647 and ONC brown dwarf systems 
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ter) by the respective authors has indicated that none fit 
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both components of the binaries simuhaneously. In the case 
of the brown dwarf EB, the models do seem to be reasonably 
consistent with both objects, but the secondary appears to 
be hotter than the primary, a very surprising result that was 
not predicted by any of the models. 

The study of low-mass stars poses a challenge for 
stellar models. Stars near to the hydrogen burning limit 
are sufficiently cool that the interior is in a partially- 
degen erate state on the main sequence (jChabrier fc Baraffd 
I1997I). and magnetic fields may play an important role 
i Mullan fc MacDonaldll200ll '). The approximations underly- 
ing the usual 'grey' atmosphere models break down, so non- 
grey model atmospheres must be determined, taking account 
of effects such as the recombination of molecules (e.g. H 2 and 
TiO) due to the low temperatures (jBaraffe et al.ll2003 ). One 
of the only ways to test these models is by comparison with 
precise observations of low-mass stars. 

It is therefore essential that a larger sample of pre-main 
sequence systems be found and characterised, to provide bet- 
ter constraints on the models. This is the primary aim of the 
Monitor project (Aigrain ct al. 2007; Hodgkin et al. 2006), 
a photometric monitoring survey we have undertaken of all 
suitable nearby, pre-main sequence open clusters and star 
forming regions, to search for detached eclipsing binary and 
transiting planet systems. This publication presents the first 
of these, a detached EB in the Orion Ne bula Cluster, whic h 
has age 1± 1 Myr and distance ~ 470 pc (|Hillenbrandlll997t ). 

Sections [2] and [3] review the photometric and spec- 
troscopic observations and analysis, respectively. Section 3] 
presents the light curve analysis and system parameters, and 
in Section[5]these are compared with the predictions of stel- 
lar models. Finally, we summarise our conclusions in SJS] 



2 PHOTOMETRY 
2.1 Survey photometry 

Eclipses in JW 380 were initially detected in our photomet- 
ric monitoring data of the ONC obtained using the 2.5 m 
Isaac Newton Telescope (INT), with the Wide Field Cam- 
era (WFC), during two 10-night observing runs, one in late- 
November 2004, and one in January 2005. This was supple- 
mented by two further 10-night runs in December 2005 and 
2006. 

A single ~ 34' x 34' field centred on the Trapezium 
region {9]j Ori) was observed for the entire time it was at 
airmass < 2.0, ~ 6 hours per night, using alternating 60 s V- 
band and 30 s i-band exposures. Due to the fast (~ 40 — 45 s) 
readout of the WFC this gave a cadence of ~ 3.5 minutes. 
A total of 1400 exposures were obtained in each passband, 
for a total time on target of ~ 80 hours. Our observations 
are sufficient to give an RMS per data point of 1% or better 
down to i ~ 17 and V ~ 18. 

Light curves were extracted for a total of 2500 ob- 
jects using our differential photometry software (|lrwin et al.l 
I2OO7I ). The V-band light curves in particular are somewhat 
affected by the presence of nebulosity in the ONC region, 
and show significantly more scatter than the i-band data. 
We therefore used the latter for detection of eclipsing binary 
systems, and the l/-band observations for confirmation. 

Due to the extensive intrinsic stellar variability seen in 
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Table 1. Photometric properties of the eclipsing binary system. 
Errors are not quoted for the INT/WFC photometry because 
these are dominated by systematic errors from the sky subtraction 
due to the extensive nebulosity in the optical. We expect that 
these are < 0.05 mag. 



ONC stars (e.g. Istassun et al.ll 19991 . iHerbst et al.ll2002h . the 
search for eclipsing systems was by necessity performed man- 
ually. The first candidate detected in this way is the subject 
of the present paper. The i-band light curve (see the up- 
per panel of Fig. [1} shows ~ 0.05 mag echpse events, with 

~ 0.03 mag peak-to-peak out of eclipse var iations. 

The object appears in the catalogue of lJones fc Walkerl 

() 19881 ) as star 380, with membership probability 0.99 de- 
rived from proper motion measurements by these authors. 
It was al so detected by the Ch andra Orion Ultradeep Project 
fCOUP: [Getman et al.ll2005l . star 468), and in the Two Mi- 
cron All-Sky Survey (2MASS), as source 2MASS J05351214- 
0531388. 

The photometric properties of the object are sum- 
marised in Table [1] V and / magnitudes were derived from 
the Monitor photometry using colour equations from the 
INT wide field survey web pagefl, and the near-IR J, 
H and K magnitudes from the 2MASS photometry us- 
ing colour equations in the 2MASS explanatory supple- 
menl[3. Spitzer/IRA C measurements are from Table 1 of 
iRebxill et al.) l|2006l l. 

The measured composite J-band magnitude and the 
models of iBaraffe et al.l l| 19981 ) imply a total system mass of 
~ 0.5 M0 if the system was a single star, or assuming that 
the primary contributes half of the system light, ~ 0.3 M0 
(the latter is a more reasonable assumption for a near equal 
mass binary), and the optical colour 1^ — / = 3.1 i mplies a 
spectr al type of ~ M5 using the intrinsic colours of iLeggetg 
1 19921 ). 

Examination of the COUP light curv^fl shows two flare 
events, but no obvious evidence for eclipses at the present 



^ 'http : //www. ast . cam. ac .uk/~wf csur/technical/photom/ colours/ | 
^ http : //www. Ipac . caltech.edu/2niass/releases/allsky/doc/sec6_4b 

Light curves are available in the source atlas on the COUP web 



site: http : //www ■ astro .psu. edu/coup/ 
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Figure 1. INT/WFC j-band light curve, plotted as a function of sample number (numbering from one for the first observation). Vertical 
lines denote the boundaries between different nights. The horizontal dashed line indicates the median light curve level. The upper panel 
shows the data, and the lower panel the data with the median light curve level subtracted from nights containing eclipses (the nights 
which were modified are shaded in the figure), to reduce the effect of the out-of-eclipse variations on the measured eclipse depths. 



time. This is under investigation, in collaboration with mem- 
bers of the COUP team. 

Spitzer/IRAC measurements for JW 380 are avail- 
able in Table 1 of iRebull et all l|2006l ), and are repro- 
duced here in Table [T] Comparing the observed [3.6] — [4.5] 
and [4.5] - [8.0] colours to those derived in NGC 2362 by 
iDahm fc Hillenbrandl 12007") indicates that these measure- 
ments may provide evidence for a weak disc excess, but this 
conclusion is somewhat ambiguous due to the obvious com- 
posite nature of the system, which will affect the observed 
colours. 

From the full INT/WFC data-set, a preliminary period 
of 2.65 days was determined using a standard box-fittin g 
least squares transit search program l|Aigrain fc Irwinl2004l ) . 



2.2 Follow-up photometry 

Initial follow-up observations were obtained during 2006 
February by J. Holtzman using the New Mexico State Uni- 
versity 1.0 m robotic telescope, in J-band, at the predicted 
times of eclipse from the INT data. 

JW 380 was also monitored using the 0.9, 1.0 and 1.3 m 
telescope^ at the Cerro Tololo Inter- American Observatory 
(CTIO). Observations were obtained from 2005 December to 
2007 January, in V and / bands, although the majority of the 
data are in /-band. Differential light curves were determined 
from PSF photometry using an algorithm for inhomogeneous 
ensemble photometr y jHonevcutt 19921) a s implemented in 
IStassun et ahl l|l999l) .[ Stassun et alj I (|2002l ) for observations 
of high-nebulosity regions such as the ONC. 

All the available data were combined to produce two 
composite light curves, one in /-band and the other in V- 
band, applying a normalisation to account for zero point off- 
sets between the various photometric systems in use, based 
on the median out-of-eclipse level. Note that the i and / 
passbands are not strictly identical: the i (SDSS) passband 
is slightly bluer than the conventional / (Cousins) passband, 
by ~ 10 per cent in the V ~ I colour. However, with the 
present data it is difficult to correct for this effect, and we 
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expect that the errors introduced by not doing so will be 
negligible given the photometric errors. 

An updated period of 2.6496 days was determined by 
applying a new double trapezoid fitting program (Aigrain et 
al., in prep). 



2.3 Out of eclipse variations 

The out of eclipse variations were found to not phase-fold 
at the same period as the eclipsing binary, and moreover 
with the full combined data-set it was not possible to find 
a consistent period for them, presumably due to changes in 
the spot coverage of the stellar surfaces causing phase and 
amplitude changes in the out of eclipse modulations, over 
shorter time-scales than the observing window. 

The INT/WFC data considered alone are concentrated 
into 4x ^ 10 night observing runs, so we attempted fitting 
of the out-of-eclipse parts of this light curve. If these are 
due to spots on the surface of one of the component stars, 
this allows the rotation period of the star to be determined, 
and if a sufficiently good model for the spot behaviour can 
be found, the modulations can be removed, to improve the 
accuracy of the eclipsing binary model fit; 

The period finding algorithm from llrwin et al.l (|2006l ) 
(based on least squares fitting of sinusoidal modulations) 
was modified to fit different phases, amplitudes and zero- 
point levels for each of the four observing runs, fitting for a 
common period, presumed to be the rotation period of the 
star. The results of this procedure are shown in Figure (2] 

Figure [2] indicates that the fit provides a good descrip- 
tion of the majority of the out-of-eclipse variations, except 
in the first five nights and the penultimate night. This in- 
dicates that a model assuming a constant stellar rotation 
period, here 4.9 days (which is not equal to the binary pe- 
riod or a multiple thereof), and with moderate evolution of 
spot coverage over one month timescales, provides a good 
description of the majority of the ou t-of-eclipse variability. 
We note also that iHerbst et al.l l|2002l ) list a rotation period 
of 4.75 days for this star, which tends to confirm our result. 

However, Figure [5] shows that the model we have pre- 
sented does not completely describe the variation in the 
nights containing eclipses, and this technique is not appli- 
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Figure 2. Fitting of the INT/WFC i-band light curve out-of-eclipse variations. The upper panels show the phase-folded Ught curve for 
each of the four observing runs, fit with the same period, but a different phase, ampUtude and zero point in each, with the fit overlaid as 
a solid line. The lower panel shows the residuals after subtracting the fit, plotted as a function of sample number. The dashed horizontal 
lines show the median light curve level. The vertical lines in the lower plot show the boundaries between nights of observations, with the 
boundaries between observing runs (corresponding to the division of the panels in the upper plot) denoted by longer black vertical lines. 



cable to the follow-up data due to the sparse time coverage 
(providing insufficient data to fit the rotational modulations 
in each individual light curve segment, before the spot con- 
figuration changes). 

Therefore, we elected to simply reduce the effect of the 
out of eclipse variability on the measured eclipse depths by 
subtracting the median out of eclipse light curve level from 
each night containing an eclipse (see the lower panel of Fig- 
ure [T]), since this technique treats the entire composite light 
curve consistently. We caution that this is not equivalent to 
modelling the spots physically, which will be necessary for a 
more accurate light curve fit to improve the derived physical 
parameter estimates. Fortunately, in most cases the effect of 
the spots on the inferred physical parameters is minor. For 
exa mple, in the ana l ysis o f the young eclipsing binary VI 174 
Ori, IStassun et all (2004 ) found that modelling the out of 
eclipse variability with spots did not alter the derived stel- 
lar parameters significantly, with changes of < 1 per cent in 
the derived masses and radii, and < 3 per cent in the ratio 
of effective temperatures. 



3 SPECTROSCOPY 

In order to confirm the eclipsing binary nature of the sys- 
tem, and to derive the orbital parameters, multiple radial 
velocity measurements are required. These were obtained 
using two instruments: the FLAMES multi-object fibre-fed 
echelle spectrograph on VLT UT2 in the optical, and the 
NOAO Phoenix spectrograph on Gemini South. 



3.1 VLT/FLAMES 

VLT/FLAMES observations were obtained during 2006 
November (programme 078.C-0841). GIRAFFE fibres were 
allocated to all of our eclipsing binary candidates in addition 
to 220 other suspected ONC members. A total of five spec- 
tra were obtained in each of two standard settings: HR15n 



covermg ~ 6450 - 6810 A with resolving power R ^ 17 000 
and HR21 covering ~ 8490 - 8980 A at E ~ 16 000. Data 
were reduced usin g the GIRAFF E Base-Line Data Reduc- 
tion Software (Ble cha et al.|[2000h . with additional custom- 
written software for sky subtraction using 13 fibres allocated 
to blank sky. The HR15n setting was strongly contaminated 
by spatially-variable emission lines from the Orion Nebula, 
so we preferred the HR21 setting for determination of radial 
velocities, and this has been used henceforth. 

Radial velocities were determined by cross-correlation 
using FXCOR in IRAF0 (|Todvl 1 19931 '). Since it is im- 
practical to observe radial velocity standard stars with 
FLAMES, we used model atmosphere s pectra from the 
R ~ 20 000 MARCS library (Gustafsso n et al.l |2003| ) to 
provide the cross-correlation templates. Cross-correlations 
for the present object were derived using a template with 
left = 3500 K, solar metallicity, and logg = 3.5, in accor- 
dance with the pre dicted surface gravity from the 1 Myr 
stellar models of iBaraffe et al.. (J_998) . Radial velocity errors 
were estimated using the method of iTonrv fc Da'^ (|l979D 
as implemented in FXCOR. 

The observed cross-correlation functions show clear 
triple-lined profiles (see Figure [S]) in a total of 7 epochs (in- 
cluding all the radial velocity data, see also i]3.2|l around 
maximum radial velocity separation, with the two outer 
components (the primary and secondary stars in the eclips- 
ing binary) exhibiting radial velocity variations. The central 
component (hereafter the tertiary - although note that this 
star is not necessarily physically associated with the binary) 
does not appear to show radial velocity variations. 

We note that strong Lil 6707.8 A absorption is clearly 
visible in the HR15n setting, which indicates youth, and 
hence membership of the ONC. The Lithium line profiles 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 3. Cross-correlation function for the first VLT/FLAMES 
spectrum (top), and a Gemini/Phoenix spectrum for comparison 
(bottom, see also i|3.2l l. showing a clear triple-lined profile in both 
cases. The central component does not appear to show signifi- 
cant radial velocity variations. A three-component Gaussian fit 
was used to determine the radial velocities, with cross-correlation 
peak heights (left to right) 0.11, 0.12, 0.23 in the VLT/FLAMES 
spectrum. 



from our 5 HR15n spectra are shown in Figure 3] Equiva- 
lent widths were estimated by fitting three-component Gaus- 
sian models to the two best-resolved epochs (the first and 
third in Fig. |3} , at the expected wavelengths from the ra- 
dial velocities for each component. The values obtained were 
0.32±0.05 A for the primary, 0.18±0.04 A for the secondary, 
and 0.12 ± 0.05 A for the tertiary (errors estimated using 
bootstrapping). Note that in order to compare to the values 
for single stars, these measured equivalent widths must be 
corrected for the relative luminosities of the stars. 

Unfortunately, we are unable to correct the Lil equiva- 
lent widths for spectral veiling at present, since we do not see 
suitable features in our spectra (e.g. Nil 6643 A, Fel 6663 A 
or Vl 6625 A as used by Palla et al.„2007i ) with which to com- 
pute the veiling, despite the majority of our spectra having 
very good signal to noise ratios (~ 40 pix~^ or 180 A~^). 
The values we have computed must therefore be considered 
to give lower limits on the measured Lithium abundance. 
Nevertheless, assuming the primary to contribute 1/2 of the 
light in the system, the implied equivalent width is ~ 0.6 A, 
which is a typical value fo r ONC stars of this mass (e.g. see 
ISiciha-Aguilar et al.|[2005[ ). 



3.2 Gemini/Phoenix 

Additional radial velocity measurements were obtained us- 
ing t he near-infrared spectrograph Phoenix (|Hinkle et al.l 
|2002| ) on Gemini South, on three nights: 2006 December 1, 2 
and 4 (programme GS-2006B-C-7). Following [Mazeh et all 
(|2002l ). we used the H6420 order sorting filter and the widest 
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Figure 4. Lithium 6707.8 A line profiles for the five 
VLT/FLAMES spectra in the HR15n setting. Arrows show the 
predicted positions of the lines from the three components (1, 2 
and 3 are primary, secondary and tertiary respectively), derived 
from the radial velocities. The primary and secondary show clear 
Lithium absorption, and the tertiary may also show it, but this 
line, if present, is not well-resolved from the primary in any of the 
spectra. 
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(4 pix ~ 0.34 arcsec) slit, with a central wavelength of 
1.558 /^m to give a wavelength range of 1.554 — 1.562 /im 
at resolving power _R ~ 35 000. Exposures were taken in AB 
pairs, nodding along the slit, to aid subtraction of the sky 
background and residual dark current features in the de- 
tector. A total of eight epochs were obtained for this object, 
with average signal to noise ratios of ~ 60 pix~^ or 220 A~^. 
Data were reduced using standard IRAF long-slit reduction 
procedures, and cross-correlations with the radial velocity 
standard star GJ 173 (Ml. 5 spectral type) in FXCOR were 
used to derive radial velocities, since this template gave the 
largest cross-correlation signal. This suggests that the pri- 
mary spectral ty pe is close to Ml. 5. Us ing the effective tem- 
perature scale of ICohen fc Kuhl (|l979l ). this corresponds to 
Toff = 3590 K. 

Despite the high resolution and signal to noise ratio of 
our data, it is difficult to see the individual line profiles in the 
Phoenix spectra. This may be partly due to the very strong 
de pendence on spectra l type in this wavelength region (e.g. 
see iBender et ani2005l ). where the atomic lines become weak 
moving to the later spectral type secondary and tertiary 
components. This also partly explains the reduction in cross- 
correlation signal for these objects (see Figure [Sj, since we 
used an earlier-type template with strong atomic lines. 

We note that the unavoidable use of different templates 
for the two instruments may introduce systematic offsets 
between the two sets of radial velocities. We have tried to 
minimise this source of error by using templates with similar 
effective temperatures. 

3.3 Spectroscopic orbit solution 

The photometry and initial VLT/FLAMES radial veloci- 
ties were found to be consistent with the photometrically- 
derived period, but the Gemini/Phoenix velocities were not. 
A consistent solution was found at twice the initial period, 
or 5.2991 days. Re-examination of the light curve indicated 
that three secondary eclipses had in fact been observed, and 
that these had depths of ~ 0.03 mag in i-band. This is sim- 
ply a result of the improved phase coverage obtained by 
including the Gemini/Phoenix radial velocity data. 

Figure [S] shows the resulting phase- folded radial veloc- 
ity curve, with period and phase zero-point fixed at the 
photometrically-determined values. Neither the light curve 
nor radial velocity curve appear to show significant devia- 
tion from the predictions for a circular orbit, with a formal 
fit giving eccentricity e = 0.004 ±0.036, so we have assumed 
zero eccentricity henceforth. 

Table [2] gives the orbital parameters derived from the 
radial velocity fit. 

The measured systemic velocity {vsys) is 26.55 ± 
0.87 km s^^. This is very close to the systemic veloc- 
ity of the ONC, of 25 ± 1.5 km s"^ (IStassun et all 1 19991 : 
ISicilia-Aguilar et al.l |2005| ) , which provides a further kine- 
matic confirmation of cluster membership. 

Since the tertiary lies approximately at the ONC sys- 
temic velocity in the cross-correlations, it seems likely that 
this star is also a member of the ONC. If this is the case, 
the m easured stellar density in the ONC (e.g. iHillenbrandl 
Il997l found ~ 1600 members with / < 17.5 over a 0.2 sq.deg 
region) indicates a probability of ^ 10"'* of a chance super- 
position of another ONC member with the binary, within 
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Figure 5. Radial velocity curve, folded on a period of 
5.2991 days, where phases 0, 1, . . . are defined to be at mid sec- 
ondary eclipse. The curves show the fit for a circular orbit, to 
the primary (solid line) and secondary (dashed line) radial ve- 
locities. The dotted line indicates the systemic radial velocity. 
Points with error bars show the radial velocity measurements 
for all three components of the system (primary: squares, sec- 
ondary: circles, tertiary: triangles), filled symbols denote radial 
velocity points from VLT/FLAMES, and open symbols points 
from Gemini/Phoenix. Velocities for the system components were 
only measured when the corresponding peaks were visible in the 
cross correlation functions, so some data points are missing for the 
secondary and tertiary components, particularly close to phases 
and 0.5. The tertiary velocities are poorly measured and are 
strongly influenced by blending of the cross-correlation profiles 
(e.g. see Figure [Sjl, so the apparent variations for this component 
are not considered significant in the present data. 



Systemic velocity 
Primary semi-amplitude 
Secondary semi-amplitude 


^sys 

K2 


26.55 ±0.87 km s"! 

32.8 ± 1.3 km s"! 

56.9 ± 2.2 km s"! 


Mass ratio (M2/M1) 


<? 


0.577 ±0.032 


Semimajor axis 


a sin i 


9.38 ± 0.27 Re 


Total mass 


M sin^ i 


0.396 ±0.019 Mq 


Primary mass 


Ml sin^ i 


0.251 ± 0.023 M0 


Secondary mass 


M2 sin'' i 


0.145 ±0.012 Mq 


Reduced 




0.60 



Table 2. System parameters derived from the radial velocity 
curve fit in Figure |5] 



0.5" (the best seeing experienced in our observing runs). It 
thus seems likely that the system is a physical triple. 



4 LIGHT CURVE ANALYSIS 

The i-band light curve was analyse d using jktebop 
llSouthworth. Maxted fc SmaUevI l2004al : ISouthworth" et al.l 
|2004b| ), a modified version of ebop (Eclipsing Binary Orbit 
Program; iPopper fc Etzcl 1981; Etzcl 1981i). These codes 
use a model in which the dis cs of the stars are ap proxi- 
mated using biaxial ellipsoids l|Nelson fc Davislll972l ). This 
approximation is only applicable in well-detached systems 
where the stars are close to spherical, as in the present case. 
Note that star spot modelling and handling of multi-band 



© 0000 RAS, MNRAS 000, 000-000 



A 0.26, 0.15 Mq pre-main- sequence eclipsing binary in the Orion Nebula Cluster 7 



light curves are not implemented in these codes, and must 
be done externally, as we have in the present work. The 
JKTEBOP code contains several enhancements, including fit- 
ting of the sum and ratio of the stellar radii, the use of 
the Levenberg- Marquardt minim ization algorithm for least- 
squares fitting (jPress et al.lll992l ). and of particular interest 
for the present discussion, Monte Carlo simulation to deter- 
mine robust error estimates for the stellar parameters, which 



has been shown to give reliable results l|Southworth et al] 
I2OO5I ). Note that we have not used the l/-band light curve 
in the present analysis, since there are no secondary eclipses 
observed, and the photometric precision reached is some- 
what poorer. The observed primary eclipses have very sim- 
ilar depth to the i-band eclipses. 

The light curve analysis for the present system is com- 
plicated by the known presence of third light. From Figure 
[S] it appears that the cross-correlation peaks for the sec- 
ondary and tertiary are of very similar height and shape, 
which suggests that the third light star has a similar flux 
and spectral type to the secondary. We therefore expect it 
to contribute ~ 1/4 of the light in the system, assuming a 
2:1:1 ratio of i-band luminosity. 

Due to the need to fit the third light contribution, we 
fixed as many parameters as possible, to improve stability. 
In 7-band, limb darkening is of much lower importance than 
for bluer passbands, so we fixed these parameters. We used 
a linear limb darkening law: 



h{i) 



= 1 - a;A (1 - 



(1) 



where fi — cos 6 is the cosine of the angle between the line- 
of-sight and the normal to the stellar surface, IxifJ^) is the 
observed surface intensity at wavelength A as a function of 
H, and x\ is the linear limb darkening coefficient. Values of 
XX of 0.7208 for the primary star (Tes = 3300 K, log g = 4.0, 
[M/H] = 0) and 0.7341 (T^b = 3100 K, logp = 4.0 
[MJH] = 0) for the secondary star were adopted from lClaretl 
l|2004h . for the SDSS i' passband, which is the closest match 
to the SDSS-like i filters used in this work. The gravity dark- 
ening exponents were also fixed, at (3 — 0.0 8, a value a ppro- 
priate for stars with convective envelopes ||LucvII1967| ). and 
the option in jktebop to calculate the reflection effect was 
used, rather than fitting for it. This is generally a reasonable 
approach for well-detached systems. 

The remaining parameters were allowed to vary. These 
are the sum of the radii, ri 4- r2 (where = Ri/a, the 
radius divided by semimajor axis, a parameter which can be 
determined from the light curve analysis alone), the radius 
ratio k = rijrx, orbital inclination i, surface brightness ratio 
J (again defined as secondary divided by primary such that 
J < 1) and fractional luminosity of the third light L3. We 
also allowed the period P and phase zero-point To (defined 
as the heliocentric Julian day of the secondary eclipses) to 
vary, to refine the existing values. Note that the quantity J 
and the luminosities quoted here are for i-band. 

Note that the ebop codes model only single-band light 
curves, and do not use model atmospheres, so it is not nec- 
essary to assume an effective temperature or surface gravity 
in the fit itself: these parameters enter only through the 
(assumed) limb darkening coefficients, and we expect this 
dependence to be weak in i-band. 

Figure|6]show the i-band phase-folded light curve, with 




Phase 

Figure 6. Phase-folded combined i-band light curve, with phase 
defined to be at mid-secondary eclipse, as Figure |3] The upper 
panel shows the light curve, with the fit overlaid (solid line), and 
the lower panel shows the residuals (data - model). The panels 
show a magnified view of the regions around primary eclipse (left) 
and secondary eclipse (right). The scatter in the secondary eclipse 
is larger than the primary since this was only observed in the 
follow-up observations and not the INT/WFC survey photometry 
(which was of better photometric precision). 



Radius sum 

Radius ratio 

Orbital inclination 

Surface brightness ratio 

Third light ratio 

Orbital Period 

Phase zeropoint (HJD) 



(Ri + R-2)la 

k 

i 

J 

L3 

P 

To 



2191+"'0035 

-+0.30 
-0.04 

ito26 degrees 
3+0.045 
-0.011 
j+0.111 
-0.058 

5.299180lg:»«li days 
2449704.45279looog4g 



0.75: 
80.42:! 
0.573:1 
0.128"* 



Primary radius 


Ri/a 


U.1248_o 0^54 


Secondary radius 


R'z/a 


0942+° °^®2 


Luminosity ratio 


L2/L1 


n Qo^ic;+0.3635 
'^■•^^"^-0.0372 


Temperature ratio 


T2/T1 


0.8700lHiS 



Table 3. System parameters derived from the light curve fit in 
Figure [6] Note that J, L3 and L2/L1 are the quantities for i- 
band. 68 per cent confidence intervals are quoted. Note that since 
the eff'ective temperatures are unknown, the temperature ratio 
was calculated assuming equal bolometric corrections for both 
components. 



the fit overlaid. The system parameters derived from the 
light curve fitting are given in Table[3l Errors were estimated 
using the Monte Carlo analysis, with 10 000 iterations, the 
results of which are also shown in Figure [T] 

The Monte Carlo algorithm uses the best flt to generate 
a synthetic light curve, injecting Gaussian noise with ampli- 
tude determined by the observational errors to produce a 
simulated light curve, which is then fit to determine a new 
set of parameters. The errors can be estimated using the 
distribution of the parameters fr om a large number of real - 
isati ons of this proces s . See a lso ISouthworth et ahl (|2004ah 
and ISouthworth et al.l ()2004lJ ) for more details. 

Since the luminosity ratio was not constrained in the 
analysis, the solution with the radii of the two stars swapped 
(A: > 1 or 7?2 > Ri) is also formally permitted by the fit. 
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0.5 1 1.5 79.5 80 80.5 81 81.5 

k = Rj / R, i / deg 



Figure 7. Confidence intervals on the light curve fit parameters from the Monte Carlo simulations. The contours show the 68.3, 95.4 and 
99.73 per cent confidence levels, corresponding to 1, 2 and 3(t standard errors. The crosses indicate the best-fitting values. The panels 
show: r2 versus ri (top left), luminosity ratio L2/L1 versus radius ratio k, fractional third light luminosity I/3 versus fc, and L3 versus 
orbital inclination i. Degeneracies between ri and r2, and L2/L\ and k (due to the existence of solutions for A: < 1 and for fc > 1 as 
discussed in the text) are clearly visible. Note that the 3o" contours are likely to be underestimated due to the limited number of Monte 
Carlo samples used. 



Primary mass 
Secondary mass 
Primary radius 
Secondary radius 



Ml 
Ma 
Ri 
R2 



(, ,^-,+0:013 

1 1CQ+0-039 p 

0.897_o;o34 R, 




M0 



•0 



Semimajor axis a 
Primary gravity log(gi) 
Secondary gravity log(g2) 



1+0.27 T, 
-0.27 -^0 
-+0.137 

-0.025 



9.52: 

3.706"* 

3.711« cms- 



Table 4. Physical parameters derived from the combined radial 
velocity and light curve fitting. Note that the errors on Mi and 
M2 , and Ri and R2 are not independent since they are calculated 
from the mass ratio and radius ratio. 



which leads to highly asymmetric errors for k, L2/L1 and 
the two fractional radii. 

Table |4] summarises the final system parameters, com- 
bining the radial velocity and light curve information. Note 
the extremely large stellar radii and low surface gravities, as 
expected for a very young system. 

We caution that the effects of the out-of-eclipse modu- 
lations are non-negligible in this system, and these have not 
been modelled at present due to the lack of suitable hght 
curve data. Therefore, the parameters derived from the hght 



curve fitting may be affected by systematic errors due to not 
properly accounting for the spot-induced photometric vari- 
ations. Improved hght curves will be required to resolve this 
problem. 



5 COMPARISON WITH STELLAR MODELS 

We are unable to determine accurate effective temperatures 
at the present time, so the only comparison which can be 
made to stellar models with the present data-set is in the 
mass-radius plane. Figure[8]shows a comparison of the obser- 
vations for the present system with the Lyon group stellar 
models. The primary and secondary are con sistent within 
the Ic r error bars with the 2 Myr model of . Baraffe et al.l 
(| 19981 ), but ~ 2(7 away from the 1 Myr model. The primary 
is also consistent within Icr with the 3 Myr model. The value 
of 2 Myr is fully consistent with t he upper end of th e canon- 
ical 1 ± 1 Myr age for the ONC (|Hmenbranl|l993). 

Figure [§] shows the position of the eclip si ng b i- 
nary relative to mo d els f ro m iBaraffe et aL 
D'Antona fc MazziteUil lll998h . ISiess et al.1 1I2OOOI) and 



Palla fc Stahlerl (Il999l). The r e sults for t he iBaraffe et al 



1998"), 'D'Ant ona fc MazzitellJ (|l998l ) and iPaUa fc StaMei 



(tl999t ) are consistent with the 2 Myr inferred age from above; 
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BCAH98 DM9B Z = 0.02, Y = 0.28, Xj=2x10' 




0.2 



Mass / 



Mass / Mg 



Figure 9. Comparison of tlie eclipsing binary system with several models of low-mass stars. The EB data are shown by the black 
points with error bars, and the lines show pre- main sequence model s for (top to bottom); 1 Myr (grey line), 2 Myr (blue), 3 Myr (red), 
5 Myr (green) and 10 Myr (magenta). Top left: iBaraffe et al] 1 1998h for solar metallicity and a = 1.0, top right: [P'Antona &: Mazzitellil 
(|l998l l for Z = 0.02, Y = 0.28, Xg = 2 x 1 0"'^, bottom left: ISiess et al.l JioOQ ) for Z = 0.02, Y = 0.277, with (moderate) convective 
overshooting, bottom right: [Palla fc Stahled l ll99£ ). The models have b een interpolated where nece ssary to provide poi nts at the required 
ages. The brown dwarf eclipsing binary system of Stassun et al.l l|2006h (using the parameters from lStassun et al.ll2007h is shown to allow 
the inferred ages to be comp ared, since this is suspe cted to also be a mem ber of the ONC, and should therefore be assigned a similar 
age to JW 380. Note that the lSiess et al.l l l2000f) and lPalla fc Stahled l ll999D models do not provide points below 0.1 Mq. 



although there are clearly visible differences between the 
results of these sets of track s (especially when comparing 
iD'Antona fc Mazzitellil Il998l to the others), the measure- 
ment errors are too large to distinguish them at present . 
Comparing the data with the models of ISiess et al.l (|2000l ) 
would indicate a slightly older age of ~ 3 Myr. However, 
there is a visible discontinuity in the latter tracks at 2 Myr, 
and they are not well-sampled over this mass range. 

Given the present observational errors, it is not yet 
meaningful to perform a more detailed comparison with stel- 
lar models. We plan to obtain improved photometry and 
radial velocity measurements to resolve this, in addition to 
spectroscopy covering suitable lines for determining effective 
temperatures and flux ratios. 



The system consists of a 0.26±0.02, 0.15±0.01 Mq eclipsing 
binary with period 5.299180 ± 0.000013 days, blended with 
a significant ('^ 20 per cent of the total luminosity) third 
light component, which appears to show no significant radial 
velocity variations over the (limited) observing window. 

Modelling the combined light and radial velocity 



curves for the system yielded stellar radii of 1.19 



-1-0.04 



and 0.90I^;o3 Rq for the primary and secondary of 
the binary system. These large rad ii are consiste n t wit h 
the 2 Myr age stella r mod els of Baraffe et al.l 
ID'Antona fc Mazzitellil (Il998l ) and iPalla fc Stahleij 



19981). 



19991 ). 



and the 3 Myr age models of ISiess et al.l l|2000l) given the 
masses of the system components. The systemic radial ve- 
locity is consistent with that of the ONC, and the detection 
of a clear Lil 6707.8 A absorption feature is also suggestive 
of youth, and hence ONC membership. 



6 CONCLUSIONS 

We have reported the detection of a new pre-main-sequence 
eclipsing binary, JW 380, which appears to be a member 
of the Orion Nebula Cluster (membership probability 0.99, 
systemic velocity within la of the cluster systemic velocity). 
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Figure 8. Mass-radius relation for low-mass stars and eclips- 
ing binaries. The present system is shown by the black points 
with error bars, and t he lines show pre-m ain-sequence NextGen 
(a = 1.0. solid lines: ISaraffe et al.lll998l) and DUSTY (dashed 
lines; IChabrier et al.ll2000h models from the Lyon group, at five 
ages (top to bottom); 1 Myr (grey line), 2 Myr (blue), 3 Myr 
(red), 10 Myr (magenta) and 1 Gyr (black). Systems shown as 
red poin ts with error bars a r e existing pre -main seque nce bina- 
ries fr o mlCovino et al.l ||200 j) , ICovino et al. Ii2004), Sta ssun et al.l 
||2004| '). iHebb et al.l ||2006|' ) and IStassun etalTteOOTl'l . The small 
grey points with er r or ba rs are a com p ilatio n of field systems 
from (Pelf ossc ct"ar (2000), ' Lane et al. | l! 200lY 'sif ;ransan et al. 



2003 ). [Lopcz-Moralcs 12004,). 



(.200, 



Pont et al. 



, . .. .Bouchy et al 
20051') and iLopez -Morales et al. I (120061) . We have opted to show 



egra 



results from the literature to produce a figure summarising the 
present empirical constraints on the PMS mass-radius relation. 
The DUSTY models are included to show the predicted behaviour 
in the brown dwarf domain. 
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